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ABSTRACT 

We quantify the active galactic nucleus (AGN) contribution to the mid-infrared (mid-IR) and the 
total infrared (IR, 8 — 1000 fim) emission in a complete volume-limited sample of 53 local luminous 
infrared galaxies (LIRGs, Ljr = 10^^ — 10^^ L©). We decompose the Spitzer Infrared Spectrograph 
(IRS) low-resolution 5 — 38 /im spectra of the LIRGs into AGN and starburst components using 
clumpy torus models and star- forming galaxy templates, respectively. We find that 50% (25/50) of 
local LIRGs have an AGN component detected with this method. There is good agreement between 
these AGN detections through mid-IR spectral decomposition and other AGN indicators, such as the 
optical spectral class, mid-IR spectral features and X-ray properties. Taking all the AGN indicators 
together, the AGN detection rate in the individual nuclei of LIRGs is ^ 62%. The derived AGN 
bolometric luminosities are in the range iboi(AGN) = 0.4 — 50 x lO'*'^ ergs""'^. The AGN bolometric 
contribution to the IR luminosities of the galaxies is generally small, with 70% of LIRGs having 
Lboi[AGN]/LiR < 0.05. Only ~ 8% of local LIRGs have a significant AGN bolometric contribution 
Lboi[AGN]/Lip. > 0.25. From the comparison of our results with literature results of ultraluminous 
infrared galaxies (-Lir — 10^^ — 10^^ Lq), we confirm that in the local universe the AGN bolometric 
contribution to the IR luminosity increases with the IR luminosity of the galaxy /system. If we add 
up the AGN bolometric luminosities we find that AGNs only account for 5%^*^ of the total IR 
luminosity produced by local LIRGs (with and without AGN detections). This proves that the bulk 
of the IR luminosity of local LIRGs is due to star formation activity. Taking the newly determined 
IR luminosity density of LIRGs in the local universe, we then estimate an AGN IR luminosity density 
of rjAGN ^ 3 X lO^L^ Mpc"^ in LIRGs. 

Subject headings: galaxies: nuclei — galaxies: Seyfert — infrared: galaxies 
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1. INTRODUCTION 

Luminous infrared (IR) galaxies (LIRGs) are defined as 
having infrared (8 — 1000 /im) luminosities in the range 
LiR = 10^^ — 10"'^^ L0. Although an active galactic nu- 
cleus (AGN) may contribute, it is believed that the bulk 
of their IR luminosity is pr oduced by dust heated b y 
intense star- forming activity (iSanders fc Mirabel II1996D . 
Using the prescription of Kennicutt (1993) the IR lumi- 
nosities of LIRGs imply star formation rates in the range 
17 - 170 Mq yr-i for a Salpeter IMF. 

The early studies of the AGN activity in LIRGs made 
use mostl y of the optical spect ral ra nge. The c lassical 
studies of IVeilleux et al. I H1995D and iKim et al. I Hl995f ) 
showed that the fraction of sources containing an AGN 
increases at higher IR luminosities, although the relative 
contributions of star formation and AGN to the bolomet- 
ric lumi nosity of the s ystem were not well determined. 
Recentlv lYuan et al. I ((2OIOD intr oduced a new approach 
for spectral classification of the IVeilleux et al 
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sample. They found that ^ 22% of LIRGs are classified 
as Seyfert or LINER, - 37% are AGN/starburst (SB) 
composites, 24% are Hll-like, and the rest are ambigu- 
ous. They also showed that a large fraction of the LIRGs 
previously classified as LINERs are AGN/SB composites. 

The AGN activity of local LIRGs has also been 
addressed using IR observations. iGoldader et al. I 
(|1997allbl ) obtained near-IR spectroscopy of a large sam- 
ple of LIRGs and foun d no new Seyfert 1 AG Ns. In the 
mid-infrared (mid-IR) iValiante et al. I (|2009D used spec- 
troscopy obtained with the Spitzer Infrared Spectrograph 
(IRS) to estimate the AGN emission at 6 fim of a sample 
of IR bright galaxies, as an ingredient for their backward 
evolution model for IR surveys. Locally they found an 
increasing contribution of the AGN emission at 6 /Ltm at 
higher IR luminosities. They did not observe this ten- 
dency at higher rcdshifts, however, when modelling the 
mid-IR and submillimeter number counts. iPetric et al. I 
([20^1 li) ■ as part of the The Great Observatories All-Sky 
LIRG Survey (GOALS, see I Armus et"!!"! 120091) . used a 
number of mid-IR spectral indicators to derive an AGN 
contribution of 12% for the GOALS LIRGs. At higher IR 
luminosities, most ultraluminous IR galaxies (ULIRGs, 
LiR — 10^^ — 10^'^ Lq) are only dominated bolometrically 
by AGN em ission at Lm > 5 x 10^^ Lq (jNardini et al. I 
I2OO8LI20T0I) . 

In this paper we decompose the Spitzer Space Tele- 
scope /IKS spectra of a complete volume- limited sample 
of LIRGs into SB and AGN components. The main goal 
is to derive the mid-IR and bolometric AGN contribu- 
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TABLE 1 

The Volume-limited Sample of Local LIRGs 



Galaxy Na 



IRAS Na 



Dist. 
(Mpc) 



log 



Rcf. 
IR 



Spo 



al Cla 



Rof. 

Class 



NGC 23 

MCG +12-02-001 
NGC 633 
ESQ 297-G012 
UGC 1845 
UGC 02982 
NGC 1614* 

CGCG 468-002-NEDOl 

CGCG 46S-002-NED02 

UGC 3351 

NGC 2369 

NGC 2388 

MCG +02-20-003 

NGC 3110 

NGC 3256 

ESQ 264-G057 

NGC 3690 

IC 694 

ESQ 320-G030 
MCG-02-33-098-W 
MCG-02-33~098-E 
IC 860 

MCG-03-34-064 
NGC 5135 
ESQ 173-G015 
IC 4280 
UGC 08739 
ESQ 221-IGOlO 
NGC 5653 
NGC 5734 
NGC 5743 
IC 4518W 
IC 4518E 
Zw 049.057 
NGC 5936 
NGC 5990 
NGC 6156 



IC 4686* 
IC 4687 
IC 4734 
NGC 6701 
NGC 6921* 
MCG +04-48-002 
NGC 7130 
IC 5179 
NGC 7469 
NGC 7591 
NGC 7679 
NGC 7769 
NGC 7770 
NGC 7771 



IRAS F00073 + 2538 
IRAS F00506+7248 
IRAS F01341-3735 

IRAS F02208+4744 
IRAS F04097+0525 
IRAS F04315-0S40 
IRAS F05054+1718 



IRAS 
IRAS 
IRAS 
IRAS 
IRAS 
IRAS 
IRAS 
IRAS 



F05414+5840 
F07160-6215 
F07256+3355 
F07329+1149 
F10015-Q614 
F10257-4339 
F10567-4310 
F11257+5850 



IRAS F11506-3851 
IRAS F12596-1529 



IRAS 
IRAS 
IRAS 
IRAS 
IRAS 
IRAS 
IRAS 
IRAS 
IRAS 
IRAS 
IRAS 



F13126+2453 

F13197-1627 

F13229-2934 

13242-5713 

F13301-2356 

F13470+3530 

F13478-4848 

F14280+3126 

F14423-2039 

F14423-2042 

F14544-4255 



IRAS F15107+0724 
IRAS F15276+1309 
IRAS F15437+0234 
IRAS 16304-6030 
IRAS F17138-1017 
IRAS 17578-0400 
IRAS F18093-5744 

IRAS F18341-5732 
IRAS F18425+6036 
IRAS 20264 + 2533 



IRAS 
IRAS 
IRAS 
IRAS 
IRAS 
IRAS 
IRAS 



F21453-3511 

F22132-3705 

F23007+0836 

F23157+0618 

23262+0314 

F23485+1952 

F23488+1949 



4478 
4722 
5128 
5183 
4601 
5355 
4778 
5268 
4951 
4433 
3196 
4078 
4908 
5014 
2790 
5141 
3057 
3098 
3038 
4714 
4775 
3859 
5009 
4074 
2948 
4843 
5000 
3084 
3513 
3998 
4121 
4720 
4518 
3858 
3977 
3756 
3253 
5161 
3931 
4948 
5105 
4623 
3896 
4327 
4199 
4837 
3363 
4840 
4907 
5162 
4158 
4128 
4277 



64.7 
68.3 
74.2 
75.0 
66.5 
77.6 
69.1 
76.3 
71.6 
64.1 
46.0 
58.9 
71.0 
72.6 
40.1 
74.4 
44.0 
44.6 
43.7 
69.1 
69.1 
55.7 
72.5 
58.8 
42.4 
70.1 
72.3 
44.4 
50.7 
57.7 
59.5 
68.2 
66.3 
55.7 
57.4 
54.2 
46.9 
74.7 
56.7 
71.6 
73.9 
66.8 
56.2 
62.5 
60.6 
70.0 
48.5 
70.0 
71.0 
74.7 
60.0 
59.6 
61.8 



11.11 
11.48 
10.64 
11.06 
11.12 
11.23 
11.67 
10.55 
10.83 
11.26 
11.13 
11.23 
11.11 
11.20 
11.66 
11.05 
11.43 
11.58 
11.22 
10.96 
10.48 
11.12 
11.00 
11.27 
11.57 
11.01 
11.05 
10.92 
10.98 
11.00 
10.84 
11.15 
10.28 
11.21 
11.00 
10.98 
11.12 
11.39 
11.31 
10.98 
11.30 
11.27 
11.03 
10.37 
10.96 
11.38 
11.18 
11.64 
11.11 
11.12 
10.90 
10.80 
11.33 



Al 
Al 
A2 
A2 
Al 



A3 
A3 



A4 
A4 



A3 
A3 



HII 

[Nov] 

Sy2 

HII 

HII 
HII 
HII 
Sy2 
LINER 
HII 



B2 
B3 
Bl 
B4 
B2 



B7 
B8 
Bl 
Bl 
B2 
B9 
BIO 
Bll 
Bll 
B8 
B2 
B2 



A2 
A2 
A3 
A3 



A3 
A3 



A2 
A2 



A2 
A3 
A3 



HU, [NoV], Sylli 
Sy2, [NcV] 



HII 
HII 



HII, [NeV] 
Sy2, [Nov] 



HII, oomposito 
Sy2, [NcV] 
[Nov] 
composito/HII 

HII 
HII 
HII 



HII, [NoV] 
Sy2:, [NeV] 

HII 
Syl, [NeV] 
oomposite, Sy2: 
Sy2/Syl, [NeV] 



B2, B5, B6, B12 



B8 

B2 



BIO 
B8 
B13 



B2 
B2 
BIO 
Bl 

B14 
314, B5, B6 



Bl, B2, B5, B6 



B2, B5, B6 



HII 



Notes. — The helioceni 
galaxies and the follow 
"no" indicates that no 
*Not observed with thi 
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tions in these systems. Section 2 presents the sample and 
the observations, while Section 3 describes the analysis 
of the data. Section 4 gives the results of the spectral de- 
composition of the IRS data into the AGN and SB com- 
ponents. Section 5 compares the AGN mid-IR detections 
with other AGN indicators including mid-IR, optical and 
X-ray indicators. Section 6 discusses the AGN detection 
rate and AGN bolometric contribution as a function of 
the IR luminosity and compares them with results for 
ULIRGs in the literature. Section 7 summarizes our con- 
clusions. Throughout this paper we assume the follow- 
ing cosmology Hq ~ 70 km s~^Mpc~^, rij\/ = 0.3 and 
0.7. 

2. THE SAMPLE AND OBSERVATIONS 
2.1. The Volume-limited Sample of Local LIRGs 



We use the volume-liniited s ample of LIRGs defined by 
lAlonso-Herrero et al. I ()2006al ). This sample was origi- 
nally drawn fr om the Revised Bright Galaxy Sam- 
ple (RBGS, ISanders et ahl [200l . which IS a com- 
plete flux-limited survey at 60 fim with flux densities 
greater than 5.24 Jy and Galacti c latitude Ibj > 5 deg. 
The distance limit imposed by lAlonso-Herrero et al. I 
(2006a) was chosen to allow for Paa observations with 
the narrow-band F190N flher of HST/NICMOS. The 
sample here has been completed to include all the IRAS 
sources in the RBGS with log(LiH/L0) > 11.05 and 
Uhci = 2750 - 5300 km s~^. For the assumed cosmol- 
ogy the distances are in the range ~ 40 — 78 Mpc, with 
a median value of 65 Mpc. 

The sample is presented in Table 1; it contains 45 IRAS 
systems. Eight IRAS systems in our sample contain 



Local Luminous Infrared Galaxies. II AGN Activity 



3 



multiple galaxies, that is, they are interacting galaxies, 
pairs of galaxies, or galaxies with companions. These 
can be readily identified in Table 1 as having the same 
IRAS name. We note, however, that MCG-03-34-063, 
whic h is part of IRAS F 13197-1627 with MCG-03-34- 
064 (jSurace et al. |[200l , is at 6394 km s"^ (from NED) 
and thus does not meet the distance criterion of our 
sample. Two galaxies, NGC 5743 and NGC 7769 (see 
Table 1), have IRS spectroscopy (see Section 2.2), but 
were not originally included in the IRAS RBGS. How- 
ever, both g alaxies are in intera ction with RBGS IRAS 
sources fsee iSurace et al. |[200l) . NGC 5743 is paired 
with NGC 5 734 and shows a diffu se Ha extension toward 
NGC 5734 (IDoDita et al. 1 12002^. NGC 7769 is part of 
the NGC 7771/NGC 7770 group, and there is evidence 
that NGC 7769 is underg oing a direct encounter with 
NGC 7771 (INordgren et alTliiggl . Additionally, these 
two galaxies not included in the RBGS have IR lumi- 
nosities similar to the rest of the galaxies (see Table 1) 
and thus we included them in our sample. The sample 
contains a total of 53 individual galaxies. 

The IR luminosities of the galaxies in the range 
8-1000 /Ltm, LiR, were calculated as defined by 
iSanders fc Mirabel I (I1996D. T he IRAS flux densities 
were t aken from lSanders et al. I ()2003[ ) and lSurace et al. I 
(|2004| ) . The latter work used image reconstruction tech- 
niques to resolve the IR emission from the individual 
galaxies of interacting systems detected by IRAS. For the 
groups of galaxies and interacting galaxies not resolved 
by IRAS we assumed that the Lir fraction of each com- 
ponent is proportional to the S'piteer/MIPS 24 ^m flux 
density (see Section 2.3) fraction. In Table 1 we list the 
IR luminosities for the individual galaxies in our sam- 
ple, rather than for the systems in the case of pairs of 
galaxies and galaxies in groups. Note that some of the 
log LiR values are below the imposed limit either because 
the galaxy is member of a pair or because of the revised 
values of the distances. 

2.2. IRS Spectroscopic Data 

All individual galaxies in our samp le, except for threqj , 
were observed by the Spitzer/mS (|Houck et al. 112001) 
instrument. We retrieved IRS spectroscopic data for 
our sample of LIRGs from the Spitzer archive. Fifteen 
/iZ^ 5 systems (16 galaxies) in our sample were observed 
in mapping mode and one in staring mode as part of 
our own guaranteed time observation (GTO) programs 
P30577 and P4047 9 (PI: G. H. Rieke) and were dis- 
cussed in detail bv lAlo nso-Her rero et al. I ()2009b[ ) and 
iPereira-Santaella et al. (2010a). Two more galaxies 
are from various programs (IC 860 from P1096 and 
NGC 7469 from P14). The res t of the sample are part 
of the GOALS legacy program (|Armus et al. Il2009| ) and 
were observed in staring mode. Observations at low- 
resolution (i? ~60-120) with the Short-Low (SL) and 
Long-Low (LL) modules and at high-resolution with the 
(i? ~ 600) Short-High (SH) and Long-High (LH) mod- 
ules were available for all the galaxies. 

For the st aring an d ma pping da ta red uction 
we followed |P ereira- Santacll a et al. I (|2010b[) and 
IPereira-Santaella et al. I (|2010al) . respectively. Briefly, 

''' IC 4686 and NGC 6921 were not observed, and the slit was 
not properly centered on the nuclear region of NGC 1614. 



we started with the basic calibrated data (BCD). Bad 
pixels were corrected using the IDL package IRSCLEAlH. 
Then we subtracted the sky emission. For the staring 
data we extracted the spectra using the standard pro- 
grams included in the Spitzer IRS Custom Extraction 
(SPICE) package provided by the Spitzer Science Center 
(SSC) and the point source calibration. For the map- 
ping observati ons we construc ted the data cubes using 
CUBISM (Smi th et al~ll2007[ ). The nuclear spectra of 
the galaxies observed in mapping mode were extracted 
using a 13'.'4x 13'.'4 aperture. Since the data cubes are 
calibrated as extended sources we applied a wavelength 
dependent aperture correction to the mapping data to 
obtain spectra comparable to those observed in staring 
mode. To calculate this aperture correction we used 
the mapping observations of standard stars (HR 7341, 
HR 6606, HR 6688, and HR 2491). 

2.3. IRAC and MIPS Imaging Data 

We used Sp itzer imaging data obtained with the IRAC 
(jFazio et al. 1 12004) and the MIPS (Rick e et al.1 120041 ) 
instruments, a s part of the GOALS legacy program 
ijArmus et al. 112009) . We retrieved the BCD from the 
Spitzer archive. The BCD processing includes correc- 
tions for the instrumental response (pixel response lin- 
earization, etc.), flagging of cosmic rays and saturated 
pixels, dark and flat fielding corrections, and fiux calibra- 
tion based on standard stars (see the IRAC and MIPS 
instrument handbooks for details). We combined the 
BCD images into mosaics using the MOsaicker and Point 
source Extractor (MOPEX) software provided by the 
SSC using the standard parameters. 

We obtained integrated photometry of our LIRGs us- 
ing an elliptical aperture chosen to encompass the emis- 
sion of each galaxy. This ellipse was calculated by fitting 
the emission of the galaxy (3tT over the background) in 
the IRAC 8/im band. In this work we only used the 
IRAC 5.8 ^m and MIPS 24 /xm images (see Section 4.2). 
For the IRAC images we applied the extended source 
aperture correction (up to 25 — 35% in the IRAC 5.8 ^m 
band) to the integrated fluxes (see the IRAC Instrument 
Handbook) to account for the diffuse scattering of the 
emission across the IRAC focal plane. 

3. ANALYSIS 

3.1. AGN+SB Spectral Decomposition of IRS Spectra 

The main goal of this work is to estimate quantitatively 
the AGN activity, both in terms of mid-IR detection and 
bolomctric contribution, in a complete sample of local 
LIRGs. In this section we flt the Spitzer /IHS SL-I-LL 
spectra with a combination of SB and AGN templates. 
Our approach to estimate the AGN contribution is simi- 
lar to other works in the literature (jNardini et al. |[2008l 
l2010f ). Basically these methods hinge on the close simi- 
larit y of the mid-IR spectra of high metallicity starbursts 
(see iBrandl et aini2006f ). and the strong differences in 
the spectral shape of AGN and SB emissions in this spec- 
tral range. This makes it possible to use templates to 
represent the SB and AGN activity (or a power-law con- 
tinuum for the latter), and thereby make a spectral de- 

^ The IRSCLEAM package is available at 

http://irsa.ipac. caltech.edu/data/SPIT ZER/i 
docs / dataanalysistools/ tools / irsclcan / 



4 Alonso-Herrero et al. 



TABLE 2 

AGN Fractional Contributions within the IRS Slit from the AGN+SB Decomposition. 



Galaxy 










NGC 23 


0.16 


0.10 


0.06 


0.02 


NGC 633 


0.15 


0.07 


0.03 


0.01 


CGCG 468-002-NEDOl 


0.54 


0.83 


0.82 


0.79 


UGC 3351 


0.25 


0.15 


0.11 


0.06 


NGC 2369 


0.26 


0.16 


0.10 


0.04 


NGC 3690 


0.81 


0.70 


0.60 


0.47 


MCG — 02-33-098-W 


0.09 


0.42 


0.33 


0.23 


MCG —03-34-064 


0.77 


0.93 


0.91 


0.85 


NGC 5135 


0.54 


0.26 


0.16 


0.09 


UGC 08739 


0.33 


0.14 


0.08 


0.03 


NGC 5653 


0.12 


0.10 


0.05 


0.02 


NGC 5743 


0.39 


0.28 


0.18 


0.09 


IC 4518W 


0.53 


0.71 


0.69 


0.60 


NGC 5936 


0.30 


0.13 


0.08 


0.04 
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NGC 6156 


0.49 


0.73 


0.65 


0.57 


IC 4687 


0.12 


0.09 


0.05 


0.02 


NGC 6701 


0.25 


0.15 


0.08 


0.04 


MCG -1-04-48-002 


0.20 


0.40 


0.41 


0.34 


NGC 7130 


0.36 


0.26 


0.18 


0.11 


NGC 7469 


0.49 


0.52 


0.43 


0.34 


NGC 7679 


0.27 


0.33 


0.24 


0.15 


NGC 7769 


0.41 


0.29 


0.19 


0.13 


NGC 7770 


0.58 


0.40 


0.27 


0.15 


NGC 7771 


0.19 


0.09 


0.05 


0.02 



composition of the observed mid-IR spectra of galaxies 
into SB and AGN components. The main difference be- 
tween our method and that of Nardini et al. is that we 
use the entire spectral range of the SL-I-LL IRS spectra, 
~ 5 — 38 iim, while they only used the 5 — 8 pm spe ctral 
range. In r e lated works both Valian te et al. I (|2009[ ) and 
iDeo et al. I (|2009[ ) derived the AGN contribution at dif- 
ferent mid-IR wavelengths of samples of local LIRGs and 
Seyferts, respectively by subtracting a scaled SB tem- 
plate from the observed IRS spectra. We will compare 
our results with theirs in Section 4.1. 

The AGN emission in the mid-IR is emission re- 
processed by warm dust in the putative dust y torus 
that surrounds the AGN (jAntonucci I I1993D . In- 
stead of representing the AGN mid-IR emiss i on as 
a power law (e.g., as done bv iNardini et al. 1 I2010H 
we used the iNenkova et al. I ()2008al |bf) clumpv dusty 
torus models together with the B ayesClumpy routine 
(|Asensio Ramos fc Ramos Almeidal l2009i see the Ap- 
pendix for more details) to fit the data. These torus mod- 
els (implemented in the CLUMPY code) are found to re- 
produce well the nuclear near and mid-IR AGN emission 
of loc a l Seyfert galaxies and PG quasars ([Niku tta ct al. l! 
20091: iMor et al. I 120091: iRamos AlmeidaeTal. i i2009i 



2011 : I Alonso-Herrero et al. 1 120111) . Moreover, clumpy 
dusty torus models in general appear to reproduce better 
the observed IR emission of local AGNs than smooth den- 
sity torus models (see e. g . lAlonso-Herrero et al. 1 120031 : 
iSchartmann et al. 1 120081: iMullanev et al. 1 120 lit ). For 
the SB compone nt we used the a verag e spectrum of lo- 
cal starbursts of 'Brandl et al. (2006*), as well as the 
LIRG templates of Riekc ct al. . (,2009i ) in the range 
log(LiR/L0) = 10.5 - 12. 



The INenkova et al. I ()2008al lbl) clumpy torus models 
are defined by six parameters describing the torus ge- 
ometry and the properties of the dusty clouds. In Ap- 
pendix Al we give a short description of the models. 
An additional parameter is the normalization needed 
to match the model to the observed IR data, which 
scales directly with the AGN bolometric luminosity 
Lboi(AGN). Our use of these models depends on the 
scaling of the torus model and the AGN bolometric lumi- 
nosity, rather than on the in dividual torus parameters. 
lAlonso-Herrero et al. I ()2011l ) recently showed that the 
values of Lboi(AGN) derived from the modeling of the 
nuclear near and mid-IR emission of local Seyfert galax- 
ies with the CLUMPY models were in good agreement 
with literature estimates. 

Details of the AGN-I-SB decomposition method are 
given in Appendix A2. Figure 1 shows a few examples 
of the results of the AGN-I-SB best fits for four different 
LIRGs with increasing AGN contributions. MCG -1-12- 
02-001 shows no evidence for the presence of an AGN, 
while the mid-IR continuum emission of IC 4518W is 
dominated by the AGN component. NGC 5135 and 
NGC 7469 are intermediate cases. The AGN is clearly 
detected in NGC 5135, and dominates the ~ 6 /im emis- 
sion within the IRS slit, but the SB component accounts 
for most of the continuum emission at A > 15 /im. The 
AGN and SB components of NGC 7469 have similar con- 
tributions within the IRS slit up to ^ 20 /im, but then 
at > 20 /im the SB component clearly takes over. The 
AGN-f-SB best fits for the rest of the sample are pre- 
sented in Figures Al and A2. For those galaxies requir- 
ing an AGN component to fit their IRS spectra. Table 2 
gives the AGN contribution within the Spitzer/IRS slit to 
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Fig. 1. — Examples of the AGN+SB best fits of the SL+LL IRS spectra for four LIRGs in our sample with increasing AGN contributions 
from top left to bottom right. The rest of the fits for the sample are presented in Figures Al and A2. For reference, the best fit model 
of MCG +12-02-001 does not require AGN emission, whereas the AGN contribution within the Spitzer/IRS slit at A > 20 ^m increases 
from NGC 5135, to NGC 7469, and IC 4518W (see Table 2). For each galaxy, the black line is the observed spectrum in the rest-frame 
wavelength, the dashed-dotted (green) line is the fitted torus model, the dashed (purple) line is the scaled SB template, and the solid red 
line is the sum of the fitted SB and AGN components. Note that for the AGN component the models only produce the continuum emission. 
Thus in the case of IC 4518W the bright AGN emission lines are not reproduced by the sum of the AGN and SB components. 



the continuum emission at four different reference wave- 
lengths relatively free of features: 6, 20, 24, and 30 /im. 

3.2. Line and Feature Measurements 

Most of the galaxies in our sample of LIRGs were 
part of the compilation of S pitzer/IRS SH and LH 
fluxes o f fine structure lines of iPereira-Santaella et alH 
()2010bD . For those not included there, the emission 
line fluxes were measured by fitting a Gaussian to the 
line plus a first or der polynomial to the local contin- 
uum (see iPereira-S antacUa ct al. 2010b, for more de- 
tails). The lines were considered as detected if the 
peak of the emission line was 3cr or move above the 
local continuum, where a is the standard deviation 
of the local continuum. The errors in the measured 
line fluxes depend on the signal-to-noise of the de- 
tected line, but are typically 10% for the bright lines 
[Ne ii] 12.81 /xm and [Ne ill] 15.56 /im. For the fainter 
lines ([O iv]25.89 ^m and [Nev]) the uncertainties are 
higher (15 — 30%), except in bright AGN as these lines 
are intense and the uncertainties are approximately 10%. 
The [O iv]25.89 /im line fluxes (or upper limits) are given 
in Table 3, together with those of the [Neil] 12.81 /im 
and [Ne III] 15.56 /im emission lines. Additionally we 
detected the [Nev] 14.32 /im line in five galaxies in 
our sa mple not included in IPereira-Santaella et al. I 
(|2010bD : CGCG 468-002-NEDOl, NGC 599 0, NGC 6156 
MCG +04-48-002, and NGC 7679 (see also lPetric et"an 



1201 It ). Table 4 lists the [Nev] 14.32 /im hue detections 
for our complete sample of LIRGs. 

To measure the apparent strength of the 9.7 /im sili- 
ca te feature frorn the SL spectra we used the definition 
of lSpoon et al. I (|2007[ ). Ssi = ln(/obs//cont) evaluated at 
9.7 /im. To estim ate the continuum emissi o n we us ed the 
same method as IPereira-Santaella et al. I (|2010al ). that 
is, we fitted a power law between the feature-free contin- 
uum pivots at 5.5 /im and 13 /im. We also measured the 
equivalent width (EW) of the 6.2 /im aromatic feature, 
also known as polycyclic aromatic feature (PAH). The 
local continuum was measured using a linear fit between 
5.75 /im and 6.7 /im, and the feature was measured be- 
tween 5.9 and 6.5 /im. Table 4 reports the values of Ssi 
and EW(6.2 /im PAH) for our sample of LIRGs. For 
objects with the silicate feature in emission Ssi is posi- 
tive, whereas it is negative for objects with the feature 
in absorption. 



RESULTS FROM THE DECOMPOSITION OF THE 
SPITZER/IRS SPECTRA 



4.1. 



AGN Dust Emission within the Spitzer/IRS Slits 

Using the AGN+SB spectral decomposition of the IRS 
low-resolution data we detected an AGN component in 
half of the LIRGs in our sample (25/50) with available 
data. Table 2 gives the AGN fractional contribution 
within the slit at 6, 20, 24, and 30 /im for each of them. 
Figure 2 shows the distributions of the AGN fractional 
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TABLE 3 
New IRS SH+LH Line Fluxes. 



Galaxy 


[O iv] 


[Ne ll] 


[Ne III] 


ESO 297-G012 


1.09 


46.68 


12.38 


UGC 02982 


1.58 


75.13 


9.07 


CGCG 468-002-NEDOl 


8.98 


9.01 


7.47 


CGCG 468-002-NED02 


< 1.07 


43.32 


5.91 


ESO 264-G057 


0.91 


46.67 


3.52 




< i.y ( 


00/1 c 
ZZ4.O0 


OV Q/1 


IC 4280 


< 0.18 


50.76 


3.76 


UGC 08739 


< 0.77 


28.85 


4.12 


ESO 221-IGOlO 


2.18 


86.74 


11.81 


NGC 5990 


11.77 


61.92 


17.86 


NGC 6156 


6.96 


35.47 


10.73 


IRAS 17578-0400 


< 0.98 


56.29 


11.46 


MCG +04-48-002 


8.42 


93.19 


20.00 


NGC 7679 


34.1 


77.29 


30.05 


NGC 7770 


< 0.8 


30.82 


7.82 



Notes. — Fluxes are in units of 10 ^^Wm . 

uum wavelengths are in good agreement. 

The fractional AGN contributions listed in Table 2 can 
also be used to get a rough estimate of the sensitivity of 
our AGN-I-SB decomposition method to the presence of 
a low-luminosity AGN. As discussed in the Appendix, 
we only considered an AGN detection when its contri- 
bution was greater than 5 — 7% at 20 [xm. Therefore 
it appears that our fitting method can only detect an 
AGN for fractional contributions C^5J^[AGN] > 0.15, 

C24Mm[AGN] > 0.05 and C|{^s^[AGN] > 0.02 - 0.04. 



Galaxy log(Cg^;;^[AGN]) 

Fig. 2. — Distributions of the AGN fractional contributions 
within the 5pit2er/IRS slit of local LIRGs at 6 /xm (lower panel) and 
24 /jm (upper panel). The AGN detections are shown as filled his- 
tograms, while the non-detections are shown as empty histograms. 
Also marked are the medians of the AGN fractional contribution 
(detections) distributions within the Spitzer/IRS slit (see Table 5). 



contributions within the IRS slit at two of the reference 
continuum wavelengths. On average the AGN fractional 
contribution within the IRS slits is relatively small, as 
can be seen from Table 2. There are only six galaxies 
with a dominant AGN contribution within the IRS slit at 
the four reference continuum wavelengths: CGCG 468- 
002-NEDOl, MCG -03-34-064, IC 4518W, NGC 6156, 
NGC 7469, and NGC 3690. 

For the sample of LIRGs the median values of the AGN 
fractional contributions tend to be similar at 6 and 20 /im 
for our sample (median C|)};?JAGN] ~ C^J}^„JAGN] = 
0.3, see Table 5), although individual galaxies show larger 
variations. For the majority of local LIRGs, the AGN 
fractional contributions decrease significantly at A > 20— 
30 jum , as also found in local Sevf ert galaxies (iDeo et al."1 
I2009D and X-ray selected AGNs IjMullanev et al. 1120111 )7 
The median value of the AGN fractional contribution 
at 24 and 30 /im for the LIRGs with AGN detections 
are C^J^s^[AGN] = 0.18 and C|{f-^s„JAGN] = 0.11, re- 
spectively. For the four Sevfert g alaxies in our sample 
in common with IDeo et aLl (|200£ 



the AGN fractional 
contributions within the IRS slit at the common contin- 



iValiante et al. I ()2009D estimated the AGN emission at 
6 /im of a sample of LIRGs by subtracting a SB template 
from the observed spectra in a much more limited spec- 
tral range, 5.6 — 6.9 /im. They e mployed the M8 2 spec- 
trum and the average spectrum of iBrandl et al. I (f2006.) . 
the latter to account for the dispersion of properties 
in star-forming galaxies. Both templates tend to have 
a lower 6 /xm cont inuum than the LIRG templates of 
IRieke et al. I ()2009D . A large fraction of the galaxies in 
our sample were included in their work. For the AGN 
detections we find a relatively good agreement, within 
10 — 50% of each other, for the AGN flux densities at 
6 /im. The largest discrepancy is for NGC 5135 for which 
we find a cl ear AGN compo nent (see Figure 1) at 6 /xm, 
whereas the lValiante et al~l (2009) upper limit is well be- 
low our estimate. Apart from our much broader spectral 
range, the differences may also arise from the different 
SB templates. We found a similar behavior from our 
fits. In those cases where w e obtained comparable fits 
with the IBrandl et al. I (j2006D average spectrum and the 
LIRG II.5OL0 template, the AGN fractional contribu- 
tion within the IRS slit at 6 /im using the latter template 
was lower, almost by a factor of t wo, than using the for - 
mer. For those galaxies for which IValiante et al. I ()2009f ) 
did not detect an AGN component, our estimated AGN 
6 /im flux densities are compatible with their upper lim- 
its. 

4.2. AGN Contribution to the Total Mid-IR Continuum 

Emission 
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TABLE 4 
AGN Indicators 



Galaxy Name 


Spectral class 


[NcV] 
Detection 


C6Mm [AGN] 


•^24^111 [AGN] 


[OIV]/[NeII] 


EW(6.2 


fim PAH) 


Ssi 


NGC 23 


composite 


No 


0.06 


0.05 


0.014 





52 


-0.36 


MCG +12-02-001 


hi: 


No 


<0.05 


<0.05 


0.016 





56 


-0.33 


NGC 633 




No 


0.08 


0.03 


0.037 





48 


-0.30 


ESO 297-G012 


hi: 


No 


<0.06 


<0.05 


0.023 





61 


-0.65 


UGC 1845 




No 


<0.07 


<0.05 


0.029 





51 


-0.94 


UGC 02982 


hi: 


No 


<0.03 


<0.04 


0.021 





57 


-0.58 


CGCG 468-002-NEDOl 


AGN 


Ygs 


0.34 


0.80 


0.997 





10 


-0.16 


CGCG 468-002-NED02 




No 


<0.08 


<0.05 


<0.025 





37 


-1.19 


UGC 3351 


Sy2 


No 


0.07 


0.07 


0.012 





42 


-1.16 


NGC 2369 




No 


0.14 


0.08 


0.015 





41 


-0.78 


NGC 2388 


hi: 


No 


<0.0S 


<0.05 


0.012 





50 


-0.51 


MCG +02-20-003 




No 




<0.05* 


<0.019 





12 


-1.47 


NGC 3110 


hi: 


No 


<0.Q4 


<0.04 


0.016 





57 


-0.33 


NGC 3256 


hi: 


No 


<0.05 


<0.04 


<0.026 





53 


-0.45 


ESO 264-G057 


hi: 


No 


<0.03 


<0.04 


0.019 





51 


-0.60 


NGC 3690 


Sy2 


No 


0.44 


0.28 


0.143 





16 


-0.60 


IC 694 


LINER 


No 


<0.04 


<0.02 


<0.064 





60 


-1.80 


ESO 320-G030 


hi: 


No 


<0.07 


<0.04 


<0.008 





51 


-0.46 


MCG — 02-33-098-W 




No 


0.05 


0.37 


<0.017 





48 


-0.02 


MCG — 02-33-098-E 


hi: 


No 


<0.0S 


<0.09 


0.039 





61 


-0.37 


IC 860 




No 


<0.05 


<0.05 


<0.173 





36 


-1.58 


MCG -03-34-064 


Syl 


Yes 


0.64 


0.85 


2.054 





01 


-0.19 


NGC 5135 


Sy2 


Yes 


0.39 


0.14 


0.658 





33 


-0.37 


ESO 173-G015 




No 


<0.Q4 


<0.04 


<0.009 





33 


-1.75 


IC 4280 


hi: 


No 


<0.02 


<0.03 


<0.004 





52 


-0.43 


UGC 08739 




No 


0.12 


0.06 


<0.027 





42 


-1.18 


ESO 221-IGOlO 


hi: 


No 


<0.02 


<0.04 


0.025 





55 


-0.32 


NGC 5653 


hi: 


No 


0.03 


0.04 


<0.018 





49 


-0.33 


NGC 5734 




No 


<0.03 


<0.04 


0.012 





42 


-0.45 


NGC 5743 


composite 


Ygs 


0.14 


0.15 


0.065 





43 


-0.57 


IC 4518W 


Sy2 


Yes 


0.34 


0.67 


2.119 





05 


-1.22 


IC 4518E 




No 


<0.04 


<0.04 


0.019 





47 


-0.56 


Zw 049.057 


composite 


No 


<0.06 


<0.05 


<0.029 





40 


-1.00 


NGC 5936 


composite 


No 


0.09 


0.07 


0.013 





50 


-0.47 


NGC 5990 


Sy2 


Yes 


0.22 


0.18 


0.190 





13 


-0.64 


NGC 6156 


AGN 


Yes 


0.05 


0.43 


0.196 





34 


0.57 


IRAS 17138 — 1017 




No 


<0.07 


<0.05 


<0.028 





54 


-0.66 


IRAS 17578-0400 




No 


<0.09 


<0.04 


<0.017 





52 


-1.10 


IC 4687 


HII 


No 


0.07 


0.05 


0.017 





63 


-0.26 


IC 4734 


HII 


No 


<0.05 


<0.05 


<0.028 





40 


-1.01 


NGC 6701 


composite 


No 


0.09 


0.06 


0.029 





50 


-0.56 


MCG +04-48-002 




Yes 


0.09 


0.41 


0.090 





50 


-0.92 


NGC 7130 


Sy2 


Yes 


0.20 


0.15 


0.186 





34 


-0.30 


IC 5179 


HII 


No 


<0.04 


<0.03 


0.017 





58 


-0.28 


NGC 7469 


Syl 


Yes 


0.32 


0.40 


0.125 





21 


-0.13 


NGC 7591 


composite 


No 


<0.04 


<0.04 


<0.010 





38 


-0.87 


NGC 7679 


Syl/Sy2 


Yes 


0.12 


0.23 


0.441 





56 


-0.05 


NGC 7769 




No 


0.06 


0.10 


0.032 





41 


-0.23 


NGC 7770 


HII 


No 


0.29 


0.27 


<0.026 





28 


-0.33 


NGC 7771 


composite 


No 


0.08 


0.04 


0.011 





45 


-0.59 



Notes. — T 

high or lov 
fractional i 



Lc spectral 
d as Seyf( 
respectiv 
sntributioi 



15 — 30% for the rest, 
silicate feature and th 



class is from Table 1. If 
rts from optical spectroj 
:ly (see Section 5.1). The 
is to the total 6 and 24 /j. 
. The EW of the 6.2 ^J.Tn 
typical uncertainties ar 



■copy , 
3 [NeV 



PAH feature 
2 0.02. *The 



issificd as HII or composite, we 
led as "AGN" or "composite" . 
n indicates if the [Nev] 14.32 fin 
!spectively. [OIV]/[NeII] is the 
I are measured in /J.m, and the 
integrated flux density is from 



list the composite cIe 
iepending on whethe 
1 high excitation line 
abserved line ratio, w 
typical uncertainties 
IRAS 25 fira converte 



ssifieation. Those galaxies i 
: the observed [Olv]25.89 fJ. 

is detected. C*^" ^ t^GN] a 
ith typical uncertai 

are 0.05 /xm. Sg^ i 
d to MIPS 24 ^m. 



ith [Nev]14.32 fim detections 
i/[NgII]12.81 fim line ratio is 
1 and C^4^^ [AGN] are the AGN 
Df 10 - 15% for bright AGN, and 
apparent strength of the 9.7 [ira 



The extent of the mid-IR emission of LIRGs shows 
no clear dependence on the IR luminosity of the sys- 
tem; some galaxies show relatively compact emission 
(~ Ikpc), while others show e xtended emission ove r 
a few kpc ([Soifer et al"! 12001 : lAlo nso-Her rero et al. 
2006bl: iDiaz-Santos et al. 1120081 : Pereira-Santaella et al. 



2010 



It is only at the high luminosity end of LIRGs 
(log(LiR/L0) > 11.80) and in ULlRGs, t hat the mid- 
IR emission starts appearing ve ry compact (jSoifer et al. I 
120001: iDi'az-Santos et al. l[2010l ). Moreover, the extent of 
the emission may depend markedly on the wavelength 
of the emission. For instance, several studies found 
that in LIRGs and star-forming galaxies the PAH emis- 
sion tends to be more extended than the mid-IR contin- 
uum (see e.g. 'Helou et al. 1 120041: lAlonso-Herrero et al. I 
[2009b; Pcrcira-Santaclla.eLalJl2010aD . 

The IRS slits do not generally cover the full extent of 
the galaxies. For the median distance of our LIRGs the 
13" slits cover approximately the central 4 kpc. In 



this section we derive the AGN contribution to the to- 
tal emission at two continuum wavelengths. To do so, 
we measured the total galaxy emission in the mid-IR 
using Spitzer imaging at two wavelengths IRAC 5.8 fim 
and MIPS 24 /im. The integrated flux densities were 
then compared with the AGN emission at these wave- 
lengths from the AGN-I-SB decomposition. Table 4 lists 
the AGN contribution to the 6 and 24 /im total emis- 
sion, C^°*„[AGN] and C^1*^„[AGN]. For the AGN non- 
detections we give the upper limits based on our esti- 
mated AGN detection limits as explained in Section 4.1. 

The median AGN contribution to the total 6 and 
24 /xm emission is 0.12 and 0.15, respectively, for those 
galaxies with an AGN mid-IR detection (see Table 5). 
Only 6% of local LIRGs have their 24 /im emission 
dominated by the AGN. Unlike what we found for the 
AGN fractional contributions within the IRS slits, for 
some LIRGs the AGN contribution to the total emis- 
sion can be higher at 24 fim than at 6 /im, especially 
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TABLE 5 

Statistics of the AGN Contribution in LIRGs with AGN 
Detections from Spectral Decomposition 



Quantity 


Mean 


Median 


16*1^ 


- 84"^ 


Within IRS Slit 




0.36 


0.33 


0.16 


- 0.54 




0.34 


0.28 


0.10 


- 0.71 




0.28 


0.18 


0.06 


- 0.65 




0.21 


0.11 


0.02 


- 0.57 


Total Emission 




0.18 


0.12 


0.07 


- 0.34 


CXmlAGN] 


0.24 


0.15 


0.05 


- 0.43 


AGN Bolometric Contribution 


Lboi[AGN]/LiR 


0.12 


0.05 


0.02 


- 0.26 



Notes. — The columns are the mean, median (50-th percentile), and 
the 16-th and 84-th percentiles. The numbers in the last column 
are the 68% confidence interval. This would be equivalent to itlo" 
of the distribution around the mean value in a normal distribution. 



for those LIRGs dominated by the AGN contribution. It 
is possible that the IRAC 5.8 /zm maps include a sig- 
nificant contributio n from the 6.2 /Ltm PAH feature (see 
iHelou et al. 1 12004 and references therein) , as well as 
some photospheric emission. Additionally, the emission 
from CLUMPY torus models is expected to be mostly 
isotropic at 24 fim, whereas that might not be the case 
at 6 i iui for obscured views (i.e., type 2) of the AGN (see 
e.g., iNenkova et "am2008b[ ). 

In summary, only in 6% (3/50) of local LIRGs is the 
mid-IR emission dominated by the AGN (C24^ni[AGN] > 
0.5). About one-quarter of the sample (12/50) have in- 
termediate AGN contributions (C^;^;;^[AGN] = 0.1 - 
0.5), while the remaining ^ 70% (35/50) have little 
(C^^J,„JAGN] ~ 0.04 - 0.1) or no AGN contribution in 
the mid-IR. 



5. COMPARISON WITH OTHER AGN DIAGNOSTICS 

5.1. Spectroscopic Activity Class 

The nuclear activity (AGN versus star formation ac- 
tivity) of local LIRGs has been the subject of numer- 



ous studies using 


optical spectroscopy (see among others. 


Veilleux et al. 1 


19951: IWu et al. 1 119981: IKewlev et al. 1 


20011; 1 Alonso-Herrero et al. 1 l2009al: lYuan et al. 1 12O1O0. 



These works made use of the classical optical spec- 
tral diagnostic diagrams ()Veilleux fc Osterbrock I [19871) 
and/or the more recent classification scheme based on 
observations of large numbers of galaxies and modeling 
(jKauffmann et al. Il2003t IKewlev et al. II2006D . Table 1 



lists the spectral activity class for the 53 galaxies in our 
sample and the relevant references. Out of the 42 indi- 
vidual LIRGs with an optical class there are: 17 pure 
Hii-like nuclei, 15 composit£0 (AGN/SB) nuclei, and 10 
AGN (Seyfert 1, Seyfert 2, LINER). Two more LIRGs 
do not have an optical class but have [Nev] 14.32 /itm 
detections and high [Oiv] 25.89 Atm/[Neii] 12. 81 /itm ra- 
tios (Table 4, and see also lPetric et al. Il2011[ ). and thus 




Galaxy log(L|-gjy-|/L0) 

Fig. 3. — Upper panel. Distribution of the [Olv]25.89^m emis- 
sion line luminosities (detections and upper limits) of our sample of 
local LIRGs. The filled histograms correspond to those LIRGs with 
an AGN detection using optical and mid-IR methods. The arrows 
mark the upper limits to the [O iv]25.89 fim luminosity in each bin. 
Lower panel. Comparison between the observed [Olv]25.89 fim 
luminosities and the IR luminosities of local LIRGs. The solid 
line is the relation expected for purely star-forming galaxies de- 
rived by Pereira-SantacUa ct al. (2010b) (see text for details), 
while the dotted lines arc the iblcr dispersion of the relation. The 
LIRG observations are color coded according to the AGN con- 
tribution to the integrated 24 ^m emission: C'24^j^[AGN] > 0.5, 

C24Mm[AGN] = 0.1 - 0.5, and C^f^^[AGN] = 0.04 - 0.1, and the 

AGN non-detections, for which C^f^^^lAGN] < 0.02 - 0.05. 



^ According to the IKewlev et al. I 1)20061 ) scheme, see also 
lYuan et al.1 l|25l0l ). 



are clearly AGNs. Finally, NGC 5743 and MCG-K04- 
48-002 were optically classified as HII, but both have 
a [Nev] detection. Both have also relatively low 
[Oiv]25.89 ^m/[Neii]12.81 ^m ratios, which are typical 
of star-forming galaxies (see Table 4 and Section 5.4). 
Therefore, although these two galaxies do contain an 
AGN, they are likely to be composite in nature. The 
rest (7 LIRG nuclei) do not have an optical classifica- 
tion. In summary, 23% of our sample are AGNs, 32% 
are Hll-like, 32% are AGN/SB composites, and the re- 
maining 13% do not have a spectral classification. 

As can be seen from Table 4, all 13 LIRGs with 
an optical classification as Seyfert nuclei and/or a 
[Nev] 14.32 /im detection have an AGN component de- 
tected through the AGN-I-SB decomposition of the IRS 
spectra. The AGN contributions at 24 /im for these vary 
from 7% for UGC 3351 to almost 90% for MCG -03-34- 
064. Among the 14 optically composite nuclei with IRS 
spectra, 8 have an AGN component detected with AGN 
contributions at 24 /xm in the range 4—28%. Three of the 
15 pure HII nuclei have an AGN component, but with a 
small AGN contribution at 24 /xm except for NGC 7770. 
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5.2. AGN Bolometric Luminosities and Comparison 
with X-ray Estimates 

The AGN bolometric luminosities derived from the 
torus model fits (Section 3.1) are between 0.4 and 
50 X lO'*'^ ergs"^, with a median value of Lboi[AGN] = 
2 X 10'*'^ ergs~^. Galaxies classified as Seyferts tend 
to show the highest values of Lboi[AGN], whereas the 
non-Seyferts have a median value of Lboi[AGN] ~ 1 x 
10'*'^ erg s~^ . The AGN bolometric contribution to the IR 
luminosity of the individual nuclei range between 1 and 
70%, with an average value of Lboi[AGN]/LiR = 12% 
(median of 5%) for those LIRGs with an AGN detection 
(see Table 5). For the AGN non-detections we estimate 
that if an AGN is present its bolometric contribution 
w ould be LhnifAGNI/ivTR < 1%. 

iPereira-Santaella et al. I (|2011| ) studied the X-ray emis- 
sion of a representative, in terms of Lir, and nuclear 
activity, sample of local LIRGs drawn from the sam- 
ple presented in this paper. Table 6 lists the eight 
Seyferts in common with their work. In this sec- 
tion we compare our estimates of the AGN bolomet- 
ric lumin osity for the Seyfert ga l axies with those ob- 
tained by IPereira-Santaella et al. I (|2011l ) from the hard 
2 — lOkeV luminosities, as well as with the values de- 
rived from the luminosity of the 6.4 keV FeKa line for 
the Compton-thick sources. The AGN bolometric lumi- 
no sities were computed u sing the bolometric corrections 
of iMarconi et al.~l (|2004[ ). All these values are listed in 
Table 6. The last column in Table 6 gives the differ- 
ences of the two estimates (taken in logarithmic units) 
of the AGN bolometric luminosity, which are on average 
0.4 dex. This is similar to the differences in bolomet- 
ric luminosity estimates found by lAlonso-Herrero et al. I 
1^01 1) for a sample of nearby Seyferts. We excluded 
NGC 7679 because it appears to be highly variable in 
X-rays (see the discussions by iDella Ceca et al. 1120011 : 
IPereira-Santaella et al~ll2011f) . 

From the upper limits of the luminosity of the hard X - 
ray FeKa emission line, IPereira-Santaella et al. I (|201lD 
concluded that in the rest of their sample, if there is an 
AGN the bolometric luminosity should be < lO*'^ ergs~^. 
This is again consistent with the typical AGN bolometric 
luminosities for the non-Seyfert LIRGs in our sample. 

5.3. The [0iv]25.89 ^m High Excitation Line 

AGNs can be readily identified in the mid-IR via 
the detection of hi gh-excitation e missi on lines (see e.g. 
Genzel et aFIHQQalSturm et al. |[200HMelendez et al. I 
200D, especially the [Nev] 14.32 /Ltm and [Nev]24.32 /xm 
lines. These are unlikely to be produced by star 
formation, but they are n ot always detected in 
relatively bright AGNs (see iWeedman et al. I 120051 : 
IPereira-Santaella et al. l[2010b( ). The detection rate of 
the [Nev] 14.32 /xm line in our complete sample of local 
LIRGs is ~ 22% (11/50), as can be seen from Table 4. 

Because the [Oiv]25.89/im emission line has a lower 
ionization potential (54.9 eV) than the [Nev] lines 
(97.1 eV), the mid-IR oxygen line is detected in a 
larger fraction of known AGNs (iPereira-Santaella et al. I 
I2010bfl . Moreover, the [Oiv]25.89 jum line appears to 
be a good indicator of the AGN power (iMelendez et al. I 
\20m iDiamond-Stanic et"an[200l iRigbv et al. II2OO90 . 
although it can also be produced in star-forming galax- 



TABLE 6 

Comparison between X-ray and AGN+SB Decomposition 
AGN Bolometric Luminosities 



Galaxy 




Lboi[AGN] 


iboi[AGN] 


Diflf. 






X-ray 


IRS fit 




NGC 3690 


3.9'' 


55. 


88. 


-0.20 


MCG -03-34-064 


15 


290. 


330. 


-0.06 


NGC 5135 


lO*' 


180. 


53. 


0.53 


IC 4518W 


2.2 


26. 


71. 


-0.44 


MCG -1-04-48-002 


6.9 


120. 


21. 


0.76 


NGC 7130 


10=* 


180. 


40. 


0.65 


NGC 7469 


17 


340. 


460. 


-0.13 


NGC 7679 


3.7/0.4*^ 


56. /4. 


75. 





Notes. — All the luminosities are in units of 10*^ ergs^^. The 2 — 
lOkeV luminosities (c orrecte d for absorption) and the Lboi[AGN] 
from X-rays are from iPereir a-Santaella et al. ( 2011). The typ- 
ical uncertainties are 5 — 10% and 25%, respectively. The for- 
mal errors for Lboi[AGN] from the IRS spectral decomposition are 
0.1 — 0.2 dex. The last column is the difference between the two 
estimates of the AGN bolometric luminosity taken in logarithmic 
units. 

''Estimated from the lumi nosity of the F eKa line. ''Variable? 
First value is from.D ella C eca et al. 1 1120011) and second value from 
IPereira-Santaella et al. 1 1 120111 ^ 




5 6 7 8 



Galaxy loglL^gj^j/Lg) 

Fig. 4. — AGN bolometric luminosity (from the AGN+SB fit 
to the IRS spectral) versus the observed [Olv]25.89^m luminosity 
for those LIRGs with an AGN detection. The galaxies are color 
coded according to the AGN contribution to the integrated 24 fim 
emission as in Fig. 3. The solid lines are the relations found for 
Seyfert Is (top, and +la dispersion dashed line) a nd Seyfert 2s 
(bott om and — Icr dispersion) in the RSA sample I IRiebv et al. I 
120091) . The error bar on the upper left corner reflects the typical 
uncertainties in the derived AGN bolometric luminosities from the 
comparison in Section 5.2, and the typical uncertainty in the [Oiv] 
fluxes (see Section 3.2). 



ies (|Lutz et al. lll998HPereira -Santa ella et al. Il2010"bl ). 

The [O iv]25.89 ^m line is detected in 70% of our sam- 
ple of local LIRGs (35/50 of the individual nuclei with 
IRS LH spectra available). In this statistic we include 
NGC 3256 and IC 694 for which the line was detected 
using IRS staring mode spectra by Berna rd-Salas et al. I 
(|2C)09.) and .Alonso-Herrcro et al. . (2009bV respectively. 
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The detection rate i n the GOALdig LIR G sample (53% 
of their sources, see iPetric et al. 11201 ll) is slightly lower 
than in our sample. By comparison the [Oiv]25.89 fim 
line is detected in ^ 26% of the ULIRGs studied by 
IVeilleux et al. I (j200 9^ , but mostly in those galaxies clas- 
sified as Seyfert and LINE R. The detection of this line 
in the iFarrah et al. 1 (|2007[ ) ULIRG sample is just under 
50%. 

In terms of the spectral class the [Oiv]25.89 /im line 
is detected in all (11) of those galaxies optically classified 
as AGN, that is, Seyferts and/or [Nev] 14.32 /im detec- 
tions. Eleven of the 17 galaxies optically classified as 
composites also show [Oiv]25.89 /im emission. Interest- 
ingly, eight LIRGs optically classified as H ii-like galaxies 
are detected in [Oiv]25.89 /im, but have no evidence of 
a hot dust component from the IRS AGN-t-SB decompo- 
sition analysis. 

Figure 3 (upper panel) shows the distribution of the 
[Oiv]25.89 /im luminosities of our sample. We also 
show in this figure as filled histograms those LIRGs 
with an AGN detection, including Seyferts, compos- 
ites, [Nev] 14.32 /im detections and LIRGs with an AGN 
component from the AGN-fSB decomposition. It is 
clear from this figure that all the LIRGs with i[oiv] ^ 
10^ L0 have an AGN detected. For the merger LIRG 
NGC 3256 we plotted here our upper limit from the 
IRS spectral mapping data, but the fiux measurement of 
iBer nard-Salas et alTT (|2009l ) from an IRS staring mode 
spectrum provides a luminosity of £[oiv] ~ 6 x 10^ L0. 
At luminosities ~ 10^ Lq < i[oiv] ^ IO^'Lq approxi- 
mately 50% of the LIRGs are found to host an AGN. 

As mentioned above the [Oiv]25.89 fim line appears 
to be a good indicator of the AGN power. For those 
LIRGs with an AGN detection we compare in Figure 4 
Lboi(AGN) from the IRS spectral decomposition and the 
observed i[oiv]- The values for LIRGs with AGN detec- 
tions are between those observed in Seyfert 1 and Seyfert 
2 ga laxies from th e Revised Shapley-Ames (RSA) sam- 
ple (jRigbv et al. 1 12009). and thus consistent with the 
[Oiv]25.89 fim emission being produced mostly by the 
AGN. 

The [Oiv]25.89 fim line can also be produced by star 
formation if the AGN contribut ion to the total luminos- 
ity of the galaxy is below 5% (jPereira-Santaella et al. I 
I2010b( ). This is indeed the case for non-AGN 
detections in our sample as we showed in Sec- 
tion 5.2. In Figure 3 (lower panel) we compare 
the observed [Oiv]25.89 fim luminosities with those 
expected from pure st a r-forma tion activity following 
iPereira-Santaella et al. I (|2010b[ ). We used the typi- 
cal [Oiv]25.89 /im/|"N eilU2.81 ^m line ra tios for high 
metallicity starbursts (|Verma et al. I [20031) and the em- 
pirical r elation between th e [Neil] 12.81 fim luminosity 
and iiR, ( Ho fc Ketol200 /) to predict the [O iv]25.89 fim 
luminosity due to star formation. As can be seen from 
Figure 3, all the galaxies with no AGN component de- 
tected through the AGN-I-SB decomposition are consis- 
tent with pure star formation. On the other hand, all the 
galaxies but one with C'24^jjj [AGN] > 0.1 are more than 
la above the relation expected for star formation. More- 
no The GOALS sample contains a total of 248 individual LIRG 
nuclei in the IRAS RBGS, and thus our sample is part of it. 



over, for all galaxies with Lroivi ^ lO^'L©, this emission 
line is mostly produced by the AGN, rather than by star 
formation. 

5.4. The [Ow]25.89 fim/[Neu]l2.81 versus EW(6.2 fim 
PAH) Diagram 

Diagnostic diagrams comparing the 

[Oiv]25.89 /im/[Neii]12.81 /im line ratio with the 
EW of one of the mid-IR PAH features are commonly 
used to assess the AG N contribution to the IR emission 
of ga laxies (see e.g. iGenzel et al. I 119981 : iDale et al~l 
120061 ). In particular, thi s diag ram has been used 

' Armus et aTl 120071 : 
and LIRGs ' 



Petric et al M2011| ). 



us giag : 

for ULIRGs '(IFarrah et al. I [2007t 
IVeiUeiix et al. 1 120091 ) 

This is because the [Neii]12.81 fim line is expected 
to be produced mostly by star formation, while the 
[Oiv]25.89 /im line tends to be more luminous in AGNs 
(jPereira-Santaella et al. Il2010bl ). In this section we will 
assess if this diagram and the AGN fraction of the total 
24 fim emission provide consistent results. 

Figure 5 presents this diagram for our sample, where 
the galaxies have been color coded according to the 
AGN contribution at 24 /im. We also show in this 
diagram a mixing curve of pure star formation and 
pure AGN emission. For the pure AGN emission we 
used the median [Oiv]25.89 //m /[Neii]12.81 fim line ra- 
tio for Seyfert galaxies 2.0 ± 1.2 (IPereira-Santaella et al. I 
I2010H) and EW(6.2 /im PAH)= 0.1 fim. For the 
pure star formation we used the average value for 
high-metallicity starbursts withou t AGN detect ions ob- 
served by ISO in the sample of iVerma et al. I ([2003) , 
[Oiv]25.89 /im/[Neii]12.81 /im=0.011. We also use the 
median value EW(6.2 fim PAH)= 0.54 ± 0.07 /im (see 
Table 4) of the Hll-like LIRGs with no AGN detection 
for the pure star formation value. 

It is clear from Figure 5 that the AGN fractions de- 
rived from this diagram are entirely consistent with the 
AGN contribution to the total 24 fim emission. For 
instance, the three galaxies with C24^jj,[AGN] > 0.5 
have [Oiv]25.89 /im/[Neii] 12.81 fim> 1 and smaU val- 
ues of the EW, which correspond to > 50% AGN frac- 
tions in the simplistic mixing curve plotted in this fig- 
ure. Most LIRGs with C^4*^„JAGN] = 0.1 - 0.5 have 
[Oiv]25.89 /im/ [Neil] 12.81 fim ratios consistent with 
AGN fractions of 5 - 50%. Those LIRGs with smaU AGN 
contributions (C24^^[AGN] < 0.1) tend to have values 
of the line ratio and the EW consistent with those ob- 
serv ed in nuclear and extra- nuclear star- forming regions 
(e.g. IDale et al. Il2006( ). On the other hand, the observed 
value of EW(6.2 fim PAH) by itself is only found to be 
a good indicator of the nuclear activity when the AGN 
contribution is dominant C2|^,jj[AGN] > 0.5. For inter- 
mediate AGN contributions there is a range of observed 
values of the EW of the 6.2 fim PAH feature, and thus 
deriving the AGN contribution is not so readily done. 

5.5. The Spoon et al. (2007) Diagram 

ISpoon et al. I (|2007l ) put forward a diagram comparing 
the EW of the 6.2 fim PAH feature and the strength of 
the 9.7 fim silicate feature to provide a general classifi- 
cation of IR galaxies. This diagram has proven useful to 
assess the presence of an AGN not only for local galaxies, 
but also for high-z IR bright galaxies (e.g.. ,Farrah et alTI 




EW (6.2/2m PAH) [fim] 

Fig. 5. — [Olv]25.89 /im/[Ne ll]12.81 fim line ratio versus EW(6.2 fim PAH) diagram for our sample of LIRGs, which are color coded 
according to the AGN contribution to the integrated 24 /im emission, as in Figure 3. The solid line is a mixing curve between pure star forma- 
tion and pure AGN emission (see text for details), whereas the dotted lines represent the Icr dispersion in the [O iv] 25.89 /im/[Ne ll] 12.81 fim 
line ratio for Seyfert galaxies and EW(6.2 fim PAH) feature for purely star-forming galaxies. The percentage numbers on the curves indicate 
the AGN contribution. 




EW (6.2/xm PAH) [/2m] 

Fig. 6. — Apparent strength of the silicate feature Ssi versus EW(6.2 fim PAH) for our LIRGs color coded as in Figure 3. If Ssi is 
positive the silicate feature is in e mission, and if Ssi is negative the feature is in absorption. The regions shown in the diagram correspond 
to those defined by ISpoon et al.' (2007) for their sample of ULIRGs. Class lA is occupied by unobscured AGNs, class IC by starburst 
galaxies (i.e., PAH-dominated galaxies), and class IB by composite (AGN/SB activity) galaxies. Classes 2A through 2C are moderately 
obscured galaxies with increasing EW(6.2 fim PAH). Classes 3A through 3C (not shown here) represent the most deeply embedded objects, 
although no ULIRGs are found in class 3C. 
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Galaxy log(L,g/LQ) 

Fig. 7. — Number of AGN detections (shaded histograms) as 
a function of the IR luminosity of the galaxy for the 53 individ- 
ual nuclei of our volume-limited sample of local LIRGs (open his- 
tograms). The galaxies have been divided in three luminosity bins. 
The light grey shaded histograms correspond to detections based on 
the IRS spectral decomposition, while the dark shaded histograms 
are LIRGs optically classified as AGNs but not detected from the 
AGN-I-SB decomposition. 



120081: IHernan-Caballero et al. II2009I) . Figure 6 shows the 
Spoon et al. diagram for our sample of LIRGs, where the 
galaxies are color- coded according to the AGN c ontri- 
bution at 24 ^m. iPereira-Santaella et al. I (|2010a|) pre- 
sented this diagram for the nuclei, integrated emission, 
and spatially resolved measurements of those LIRGs in 
our sample observed in spectral mapping mode. 

On the Spoon et al. diagram a large fraction (40%) 
of the LIRGs appear mostly concentrated in the re- 
gion occupied by starburst galaxies or class IC. Approx- 
imately 25% are in class IB, which is where compos- 
ite (AGN/SB activity) galaxies are located. Only one 
LIRG is in the region of unobscured AGNs (Seyfert Is 
and quasars) or class lA, which typically tend t o show 
low values of the EW of the PAH features (Roche et al. I 
1991; Rig opoulou et al. 1 119991: iBuchanan et al. 1 120061: 
Tommasin et al. I I2008D . The rest of the sample, ap- 



proximately 30%, lie in those regions occupied by mod- 
erately obscured nuclei (—2 < Ssi < —0.8), classes 2B 
and 2C. Only one LIRG is in region 2 A. There are no 
deeply embedded {Ssi < —2) nuclei among the LIRGs 
in our sample. On this diagnostic diagram ULIRGs are 
mostly located along two branches. The first one is a 
horizontal branch that goes from class lA to class IC, 
and the second is a diagonal going from class 3A, which 
are d eeply obscured nuc lei (not shown here), to class IC 
(see iSpoon et al. 1120071 for more details). 

The Spoon et al. diagnostic diagram for LIRGs only 
provides a relatively good separation between AGN- 
dominated LIRGs, LIRGs with intermediate AGN con- 
tributions and LIRGs with very little or no AGN contri- 
bution. Only one of the galaxies dominated by the AGN 
emission (C24^,„[AGN] > 0.5) appears in the lA class, 
whereas the other two are in the composite region and 
in the 2 A class. Most of those with intermediate AGN 
contributions (Cj^i^m [AGN] = 0.1 - 0.5) are in the com- 
posite region, as expected. However, two of them appear 



in region IC, which is that occupied by PAH-dominated 
galaxies. Those with little {C^f^^[AGl^] = 0.1-0.04) or 
no AGN contribution are mostly either in the IC region 
(starburst region), or near the border between class 2B 
and class 2C, and therefore are not easily identified as 
having an AGN from this diagram. 

It is interesting to note here that the iNenkova et a"l~l 
(|2008ai rbl) clumpy torus models never produce a very 
deep silicate feature (i.e., the models have Ssi ^ — !)■ 
This means that for deeply embedded g alaxies the AGN 
temp lates used here (see discussion bv ILevenson et al~l 
120071 ) may not be appropriate and an AGN component 
might be missed. For instance two galaxies with rel- 
atively deep silicate features, IC 694 and MCG-f02-20- 
003, are classified as LINER and composite, respectively, 
but an AGN component has not been detected through 
our spectral decomposition fit. However, it is impor- 
tant to note here that deeply embedded sources are rare 
among local LIRGs. 

6. DISCUSSION 
6.1. AGN Detection Rate in Local LIRGs 

In this section we summarize the results of Sections 4 
and 5 in terms of the AGN detection rate in our com- 
plete volume-limited sample of local LIRGs and inves- 
tigate any possible dependence with the IR luminosity. 
We also compare the ACN detection rate in LIRGs with 
other types of galaxies. 

Out of the 50 galaxies with IRS spectroscopy we de- 
tected an AGN component in 25 using the AGN-I-SB 
spectral decomposition (see Table 2, and Figures 1 and 
Al). An additional eight galaxies are classified spectro- 
scopically as composite or LINER. Of these, one did not 
have IRS spectroscopy (the IRS spectrum of NGC 1614 
was not well centered) and for the remaining seven the 
AGN was not detected based on our spectral decomposi- 
tion (UGC 1845, MCG -t-02-20-003, IRAS 17138-1017, 
IC 694, NGC 5734, NGC 7591, Zw049.057). This brings 
the combined optical+IR AGN detection rate in the in- 
dividual nuclei of our sample of LIRGs to ^ 62%. In 
terms of the 45 IRAS systems, an AGN is found in 32, 
that is, ^ 70%. This AGN detection rate is in good 
agreement with tha t derived from optical emission lines 
(I Yuan et al. IIMol) and X -ra y emi ssion (jRisaliti et al. I 
l2000 HPereira-Santaella et al. II20T1I) . 

Figure 7 shows the number of AGN detections for the 
individual nuclei of our sample of LIRGs in three IR lumi- 
nosity bins. As can be seen from this figure, the AGN de- 
tection rate stays approximately constant at a 60—65% 
rate for the three luminosity bins log(LiR/LQ) < 11, 
\0 R(Lm/Lr:,) = 11-11.3 and log(LiR/L0) = 11.3-11.7. 

IPetric et al. I (j2011 ) have conducted a study similar to 
ours using the GOALS sample of 248 LIRGs. They iden- 
tified AGN through [Nev] 14.32 /^m detections in 18% of 
their sample, virtually identical to the 22% detection rate 
for us. They also detected [O iv]25.89 ^m in 53% of their 
galaxies, whereas we detected it in 70% (see Section 5.3). 
They do not count the [O iv] detections as AGN, whereas 
we find that many of these galaxies do indeed harbor ac- 
tive nuclei (from the spect ral decomposi t ion an d / or opti- 
cal class, see Section 5.3). IPetric et al. I (|2011[ ) also used 
the EW of the 6.2 /im aromatic feature and the 5 to 15 /im 
continuum flux ratios as AGN indicators, and concluded 
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Fig. 8. — Left panel: Comparison between the AGN fraction at 24 /xm and the ratio of the AGN bolometric luminosity to the IR 
luminosity of the galaxy for our sample of local LIRGs. The symbols are size (and color) coded according to the IR luminosity of the 
galaxy. Right panel: A similar comparison but at 6 /im. 



that ~ 10% of their sample are dominated by AGN in 
the mid-IR. We find that ^ 6% (Table 4) of our sample 
have a mid-IR AGN contribution greater than 50% of 
the total luminosity at 24 /im. The difference may arise 
from slightly different methodologies, but both studies 
agree that the AGN plays a substantial role in the mid- 
IR energetics only in a relatively small minority of these 
galaxies. However, we also find that 50% of the galax- 
ies in our sample contain AGNs, based on our mid-IR 
spectral decomposition method, and ~ 60% if we com- 
bine optical and mid-IR indicators. Our ab ility to find 
many more AGNs than iPetric et al~l ()201lD did in the 
larger but similar sample suggests that our spectral de- 
composition method is a powerful way to identify subtle 
AGN features. As a result, we can show that AGNs ac- 
company star formation activity in a large proportion of 
local LIRGs, although in most cases the AGNs are not 
energetically important. 

The AGN detection rate in local LIRGs is very similar 
to that of local ULIRGs. The latte r is found to be 70% on 
average using m id-IR diagnostics (jlmanishi et al. I [20071 : 
iNardini et al. | [2010). In lower luminosity galaxies AGNs 
are detected in much smaller fracti ons. For instance, 
in the optically-selected RSA sample iMaiolino fc Riekel 
(|1995D found a 5 — 16% AGN incidence, based on optical 
indicators. The AGN fraction in nearby {d < 15 Mpc) 
moderately IR luminous galaxies (Lir > 10^ L©) is 
slightly higher, 27%, when in cluding [Nev] 14.32 ij,m 
detect ions, as demonstrated by iGoulding fc Alexander I 

6.2. AGN Bolometric Contribution in Local LIRGs 

Although the combined optical-fIR AGN detection 
rate in local LIRGs is relatively high, ~ 62%, and only 
slightly less than in local ULIRGs, the important quan- 
tity is the AGN bolometric contribution to the IR lu- 
minosity of the galaxies. Using the AGN bolometric lu- 
minosities derived from the spectral decomposition (Sec- 
tion 5.2), we can derive the AGN bolometric contribu- 



tion to the total IR luminosity of local LIRGs (both with 
and without AGN detections). We find that AGNs con- 
tribute 5tl% of the IR luminosit}E3. The upper and 
lower limits take into account the uncertainties in the 
derived AGN bolometric luminosities. For the 9 galaxies 
with an AGN optical classification but no AGN emission 
detected from the AGN-I-SB decomposition, we assumed 
the typical AGN bolometric luminosity of the composite 
sources (see Section 5.2). This estimate proves that the 
bulk of the IR luminosity of local LIRGs is due to star 
formation activity. 

Our estimate of the AGN contribution to the total IR 
luminosity of LIRGs is in excellent agreement with that 
estimated from the X-ray propertie s of the rep resenta- 
tive sample studied by Pereira- Santa ella et al. I (|2011). 
IPetric et ai~l ()2011[ ) estimated that AGN supply 12% 
of the total energy emitted by the GOALS LIRGs, which 
is compatible with our upper limit of the AGN contribu- 
tion. 

Taking the newly determined IR luminosity density 
produced by LIRGs in the local universe of Oir = 
6 X lO^LoMpc"^ (jGoto et al. II2011I ). we estimate that 
the AGN IR luminosity density in the LIRG luminos- 
ity range is ^^^^ = itl x IO^LqMpc"^. Our estimate 
of the AGN IR lu minosity density is about seven times 
lower than that of iGoto et al. I (|2011[ ). probably because 
these authors assumed that in those LIRGs classified as 
active galaxies and especially those classified as compos- 
ite the IR luminosity is entirely due to the AGN. 

Finally, we can also compare Cg°*j„[AGN] and 
C^l^-al-^GN] with the AGN bolometric contribution to 

the IR luminosity of the galaxies (Figure 8). The re- 
lation between the AGN contribution in the mid-IR 
and the AGN bolometric contribution is slightly bet- 

Strictly speaking, this is just an upper limit to the AGN 
contribution because we are assuming that in LIRGs Lijoi{AGN) ~ 
Lir(AGN). 
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ter at 24 fim than at 6 fim. As explained above, the 
IRAC 5.8 fj,m images of some LIRGs may include a siz- 
able contribution from the 6.2 fim PAH feature, and/or 
the 6 fim emission of type 2 AGNs might not be fully 
isotropic. It can also be seen from Figure 8 that, al- 
though C^f^„JAGN] and Lboi[AGN]/LiR do not follow 
a 1:1 relation, the relation is nearly linear. On aver- 
age C*2J„„[AGN] = (2.9 ± 1.6) X Lboi[AGN]/LiR. We 
conclude that estimating the AGN fraction at 24 /im is 
a good proxy for the bolometric AGN contribution in 
LIRGs. Moreover, because the fractional AGN contribu- 
tion is on average higher in the mid-IR than bolomet- 
rically, the mid-IR spectral range is very appropriate to 
study the fractional AGN/SB contribution. Wu ct al. 
(j2011| ) reached a similar conclusion for a sample of IR 
luminous galaxies at z 0.3. 

6.3. Comparison with Local ULIRGs and high-z IR 
bright galaxies 

We summarize the AGN bolometric contribution for 
our sample of local LIRGs (with and without AGN de- 
tections) in Figure 9. We also compare in this figure our 
es timates for local LIIl Gs with those for local ULIRGs 
of lNardini et"!!"! (|20Tol) . This comparison is meaningful 
because both our work and that of Nardini et al. quantify 
the AGN bolometric luminosities. In this figure we sepa- 
rated the AGN bolometric contributions iboi [AGN] / Lir 
in three ranges: < 0.05, 0.05 - 0.25, > 0.25. In ap- 
proximately one-third of local LIRGs the AGN bolomet- 
ric contribution is mild, Lboi [AGN] /Ljr > 0.05. Only 
~ 8% of local LIRGs have a significant AGN contribu- 
tion, Lboi[AGN]/LiR > 0.25. 

It is clear from Figure 9 that the fraction of galax- 
ies with Lboi[AGN]/LiR > 0.05 increases at higher Ljr, 
going from 30% in local LIRGs to ~ 50% to ULIRGs 
with LiR < 3 X 10^^ L0. However, it is only at _Lir > 
5 X 10^^ L0 that the AGNs start dominating bolomet- 
rically the IR luminosity in a large fraction of local 
ULIRGs; 40% of ULIRGs in the high luminosity bin in 



Figure 9 have Lbni [ AGN] /Ltr > 0.60 (see Figure 9 of 
INardini et an[2010h . 

In summary, in the local universe there is an increasing 
bolometric contribution from AGNs at higher IR lumi- 
nosities, goi ng from 5% in LIRG s to an average of 27% 
in ULIRGs (INardini et al. l[2010l) an d 35 - 40% for the 
IJy ULIRG sample (|Veilleux et al. 1 120091. 

While the trend for an increased AGN bolometric dom- 
inance at high IR luminosities is well established locally, 
at high-z it is still a matter of debate. For instance, 
iValiante et al. I ()2009l ) needed to introduce an evolution 
in the AGN contribution to model the submillimeter to 
mid-IR number counts and redshift distribution of high- 
z IR galaxies. In particular they found that at a given 
IR luminosity high-z IR galaxies needed a smaller AGN 
contribution t han lo cally. This is in line with results of 
iFadda et al. 1 (l20Toh using deep mid-IR spectroscopy of 
z ~ 1 LIRGs and z '--^ 2 ULIRGs. In both populations 
the fraction of AGN dominated sources is si nall, being 
~ 12% in ULIRGs and ~ 5% in LIRGs. Also lWu et all 
(pOnl ) found no significant change in the overall star for- 
mation contribution to Lir from LIRGs to ULIRGs at 
z ^ 0.3. However, other works found evidence for higher 
AGN contributions for the most IR luminous ga.laxies 
at high-z. Indeed, iMenendez-Delmestre et al. I ()2009l ) 
found that, although submillimeter galaxies are domi- 
nated by intense star formation, the average AGN bolo- 
metric contribution could be as high 30%. Moreover, 
iRuiopakarn et al.~l (|201lD found a strict upper limit of 
~ 5 X 10^^ — 10^"^ L0 for the IR luminosity to be due to 
star formation in galaxies at all redshifts. Above that 
limit, AGN are expected to provide the power. 

7. CONCLUSIONS 

We have decomposed the Spitzer/IKS ~ 5 — 38 ^m 
spectra of a complete volume-limited (d ~ 40 — 78 Mpc) 
sample of local LIRGs into AGN and SB components. 
The main goal of this work is to quantify the AGN con- 
tribution to the mid-IR and IR emission of these sys- 
tems. For the SB component we used starburst and 
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LIRG te mplates, whereas t he AGN component was fitted 
with the iNenkova et al. I ()2008airbl ) clumpy torus mod- 
els. We have also compared our mid-IR AGN detections 
with this method with other AGN indicators, including 
mid-IR spectral features, optical classification and X-ray 
properties. Our main conclusions can be summarized as 
follows, 

1. Using the AGN-I-SB spectral decomposition we de- 
tected an AGN component in 25 out of the 50 in- 
dividual LIRG nuclei with IRS spectroscopy. For 
an additional nine galaxies optically classified as 
composite (AGN/SB) we did not detect an AGN 
with our method. The combined optical and mid- 
IR AGN detection rate is ~ 62% (33/53) for the 
individual nuclei of local LIRGs. This is in good 
agreement with the AGN detection rate obtained 
from optical spectroscopy, if composite objects do 
contain an AGN. 

2. The AGN contribution to the mid-IR continuum 
emission within the IRS slits is small, and decreases 
toward longer wavelengths. The median values 
range from C^)}^ [^qjs^j ^ q.SO, and C^J^,„[AGN] = 

0.18 to Cll^s jAGN] = 0.11. 

3. We used IRAC 5.8 fim and MIPS 24 /xm im- 
ages to derive the AGN contribution to the total 
emission at these mid-infrared wavelengths. The 
median AGN contributions for those LIRGs with 
an AGN detection are C^°*„[AGN] = 0.12 and 
C24lm[^GN] = 0.15, with only - 6% (3/50) of 
local LIRGs having a dominant AGN contribution 
at 24 fim (C^f^„[AGN] > 0.5). 

4. We detected the [Oiv]25.89 /im high excitation 
emission line in 70% (35/50) of the individual LIRG 
nuclei. All the galaxies in our sample with L[oiv] ^ 
10^ Lq contain an AGN and these luminosities are 
consistent with those expected from their derived 
AGN bolometric luminosities. On the other hand, 
the [Oiv]25.89 /im luminosities of those galaxies 
without an AGN detection can be explained as pro- 
duced by star formation activity. 

5. When compared with other mid-IR spectral AGN 
diagnostics, we found that our 24 /im AGN frac- 
tional components are consistent with those de- 
rived from the [OIV]/[NeII] versus EW(6.2 /im 
PAH) diagram. Using the Spoon et al. diagram 
we can identify the presence of an AGN when 
Ct2U[AGN]>0.1. 

6. From the scaling of the fitted torus models to 
the AGN component of the IRS spectra we de- 
rived AGN bolometric luminosities in the range 



Lboi(AGN) = 0.4-50 x lO'^^gj-gg-i ^j^h a median 
value of Lboi( AGN) = 2 x lO'^^ g^gs-i. These are in 
a fairly good agreement with those estimated from 
hard X-ray measurements after applying a bolo- 
metric correction. 

7. One-third of local LIRGs have Lboi[AGN]/LiR, > 
0.05, with only ~ 8% having a significant AGN 
contribution Lboi[AGN]/iiR > 0.25. This is in 
contrast with the ~ 20% of local ULIRGs with 
LiR < 3 X 10^2 L© and the ~ 60% of ULIRGs with 
LiR > 3 X 10^2 L© having iboi[AGN]/LiR > 0.25, 
respectively. 

8. Adding up the AGN bolometric luminosities in our 
sample we find that AGNs are only responsible for 
5^3% of the total IR luminosity produced by local 
LIRGs. This translates into an AGN IR luminosity 
density of rjfjf^ = 3t^ x IO^Lq Mpc"^ in local 
LIRGs. Our results prove that the bulk of the IR 
luminosity of local LIRGs is due to star formation 
activity. 

In summary, mid-IR spectral decomposition is a pow- 
erful tool to estimate the AGN contribution to both the 
mid-IR emission and the total emission of IR-selected 
sources. This is because not only are we able to identify 
all Seyfert-like AGN in local LIRGs but also because this 
technique is powerful enough to identify subtle AGN fea- 
tures (i.e., low luminosity AGN) that might be missed by 
other mid-IR methods. As a result, we have shown that 
AGNs accompany star formation activity in the large ma- 
jority of local LIRGs, although in most cases the AGNs 
are not energetically important. 
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APPENDIX 

Al. THE CLUMPY TORUS MODELS 

We used the CLUMPY models ("Nenk ova et al. |[2008allH ) to represent the AGN continuum emission in the mid-IR. 
These models are described by six parameters that deal with the torus geometry and cloud properties. These are: (1) 
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the radial thickness of the torus Y = Ro/ R^, where i?o and i?d are the outer and inner radii of the torus, respectively. 
The inner radius of the torus in these models is set by the dust sublimation temperature, which is assumed to be 
Tsub ~ 1500 K. (2) The angular distribution of the clouds, which is assumed to have a smooth boundary, is described 
as a Gaussian with a width parameter CTtorus- (3) The radial distribution of the clouds, which is described as a declining 
power law with index q (cx r""?). (4) The mean number of clouds along a radial equatorial ray Nq. (5) The optical 
depth of the clouds Ty. (6) The viewing angle to the torus, i. In the models the radiative transfer equations are solved 
for each clump and thus the solutions depend mainly on the location of each clump within the torus, its optical depth, 
and the chosen dust composition. For t he fits done in this work we adopted a dust extinction profile corresponding to 
a standard cold oxygen-rich ISM dust (jOssenkopf et al. II1992D . The total torus emission is calculated by integrating 
the source function of the total number of clumps convolved with the radiation propagation probability along the torus 
(|Ncnkova ct al. 2002). For unobscured views of the AGN, it is also possible to include its contribution to the resulting 
IR emission. The AGN co ntinuum emission in these models is characterized with a piecewise power law distribution 
(see lNenkova et al~ll2008al for details). In addition to the the six torus model parameters, there is an extra parameter 
to account for the vertical displacement needed to match the fluxes of a given model to the observations. Th is vertical 
shift, which we allow to vary freely, scales with the AGN bolometric luminosity fsee lNenkova et al. Il2008bl ). 

Given the large number of CLUMPY models (currently more than 10^) and the inherent degeneracy in the torus 
model parameters, lAsensio Ramos fc Ramos Almeidal ()2009D took a Bayesian approach to fit models to the data and 
to derive meaningful confidence levels for the fitted parameters. To this end, they developed a tool called BayesClumpy 
that allows fitting both photometric points and/or mid-IR spectra. For the Bayesia n inference of BayesClumpy we used 
an int erpolated version of the CLUMPY dnsty torus models. We refer the reader to lAsensio Ramos &: Ramos Almeidal 
(|2OO90 for details on the interpolation methods and algorithms used by BayesClumpy to perform the torus model fits 
to the data. 

A2. ACN+SB SPECTRAL DECOMPOSITION FITS 

We used an iterative method to perform the AGN+SB decomposition. In the first step we st arted by using onl y two 
AGN templates given the large number of CLUMPY torus models and their degeneracy (see iRamos A lmeida et al. I 
IIM I2011I : E sensio Ramos fc Ramos Almeidal I2009L for detailed discussions on this issue). We used the best 
fit CLUMPY torus inodel o f the average spectral energy distributions of Seyfert Is and Seyfert 2s inferred by 
IRamos Almeida et al. I (|2011| ). We obtained the best fit to the IRS spectra for each LIRG in our sample with this 
initial combination of AGN+SB templates. We next subtracted the SB component from the original IRS spectrum 
to obtain the AGN-only mid-IR sp ectrum. This AGN-only spectrum was then fitted with the BayesClumpy routine 
(jAsensio Ramos fc Ramos Almeida 2009) to infer the best iit torus model and the corresponding scaling for each 
galaxy. We finally refitted the observed IRS spectrum of the LIRG using its corresponding best-fit torus model, allow- 
ing for rescaling, and the appropriate SB template that minimized . The fits were done in the 6 — 30 /J,m spectral 
range to avoid the edges of the spectra and the slightly decreased signal-to-noise ratio of some of the spectra observed 
in mapping mode at A > 30 /xm. We note, however, that the fits tend to be good also at A > 30 ^m. We considered 
an AGN detection when its contribution was greater than 5-7% at 20 /im. Below this limit we were not able to fit 
the AGN-only spectra with BayesClumpy. This is because after subtracting the SB template in these cases a large 
proportion of the data points were close to zero or even negative. For each fit the formal errors of the AGN fractional 
contribution within the IRS slit were calculated using the observational errors of the IRS spectra and the covariance 
matrix. The typical errors of the AGN fractional contribution within the IRS slit at 6 fxm are 0.01 — 0.03. If the models 
fit the data, as is the case for most galaxies, these can be taken as the typical uncertainties of the AGN fractional 
contributions. This is because the variation of the fitted parameters is compatible with the data observational errors. 
The uncertainties associated with the use of different SB templates are discussed briefly in Section 4.1. 

We present the best AGN-I-SB fits for all the LIRGs in our sample in Figures 1, Al, and A2. As can be seen from 
Figure A2 a few LIRGs deemed not to have a mid-IR AGN detection have best fit requiring a small AGN contribution. 
These are NGC 5734, which is optically classified as composite, and ESQ 221— GOlO and IC 4734, which are classified 
as HII (see Table 1). The AGN fraction at 20 /im for all of them was, however, below the imposed limit. In these 
cases we were not able to fit the AGN-only emission with the torus models, and the AGN template corresponds to the 
Seyfert 2 template used in the first step of the iteration. 
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Figure A2. AGN+SB decomposition fits for those LIRGs without a mid-IR AGN detection. Lines are as in Figure 1. 
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